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Microcanonical Thermodynamics [1] allows the application of Statistical Mechanics on one hand to 
closed finite and even small systems and on the other to the largest, self-gravitating ones. However, 
one has to reconsider the fundamental principles of Statistical Mechanics especially its key quantity, 
entropy. Whereas in conventional Thermostatistics the homogeneity and extensivity of the system 
and the concavity of its entropy S{E) are central conditions, these fail for the systems considered 
here. E.g. at phase separation the entropy S{E) is necessarily convex to make e'^'^)"^''^ bimodal 
in E (the two coexisting phases). This is so even for normal macroscopic systems with short-range 
coupling. As inhomogeneities and surface effects in particular cannot be scaled away, one has to 
be careful with the standard arguments of splitting a system into two or bringing two systems into 
thermal contact with energy or particle exchange [2]: Not only the volume part of the entropy 
must be considered. The addition of any other external constraint like a dividing surface reduces 
the entropy. As will be shown here, when removing such constraints in regions of a negative heat 
capacity, the system may even relax under a flow of heat (energy) against the temperature slope. 
Thus Clausius formulation of the Second Law: "Heat always flows from hot to cold" can be violated. 
Temperature is not a necessary or fundamental control parameter of Thermostatistics. However, 
the Second Law is still satisfied and the total Boltzmann-entropy keeps constant or rises when 
constraints are removed. In the final sections of this paper the general microscopic mechanism 
leading to condensation and to the convexity of the microcanonical entropy S{E) at phase separation 
is sketched. Also the microscopic conditions for the existence or non-existence of a critical end- 
point of the phase-separation are discussed. This is explained for the liquid-gas and the solid-liquid 
transition. 



I. WHAT IS ENTROPY? 

Entropy (S) is the key quantity of all Thermostatis- 
tics and Thermodynamics. Therefore, its proper under- 
standing is essential. In the literature this is sometimes 
clouded by the frequent use of the thermodynamic limit 
and/or the use of the Boltzmann-Gibbs canonical statis- 
tics. In what follows, I will try to define it with minimum 
bias. 

Boltzmann's epitaph reads 



appropriate for a finite quantum-mechanical system: 



S=k*lnW 



(1) 



with 



r d^^ P d^^ q - - 

WiE^N^V)^ J ^,(2,^)3/ ^o^(i?-g{g,P}) (2) 

in semi-classical approximation. E is the total energy, N 
is the number of particles and V the volume. Or, more 
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W{E,N,V) = Tr[PE] 



(3) 



E 



all cigcnstatcs n of H with given N,V, 
and _E - Co/2 < B„ < -E + eo/2 



and eo « the macroscopic energy resolution. This is still 
up to day the deepest, most fundamental, and most sim- 
ple definition of entropy. There is no need of the ther- 
modynamic limit, no need of concavity, extensivity and 
homogeneity. In its semi-classical approximation, eq.(2), 
W{E, N,V,- ■ ■) simply measures the area of the sub- 
manifold of points in the 6 A^-dimensional phase-space (T- 
space) with prescribed energy E, particle number N, vol- 
ume V, and some other time invariant constraints which 
are here suppressed for simplicity. Because it was Planck 
who coined it in this mathematical form, I will call it the 
Boltzmann-Planck principle. 

There are various reviews on the mathematical foun- 
dation of Statistical Mechanics e.g. the detailed and in- 
structive article by Alfred Wehrl [2]. He shows how the 
Boltzmann-Planck formula (1), with (3) can be gener- 
alized to the famous definition of entropy in Quantum 
Mechanics by von Neumann [3] 

S ^ -Tr[p\n{p)], (4) 

addressing general (non projector like) densities p. 
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Wchrl discusses the conventional canonical Boltzmann- 
Gibbs statistics where all constraints are only fixed to 
their mean, allowing for free fluctuations. He points to 
the many serious complications of this definition. How- 
ever, in the case of conserved observables we know more 
than their mean. We know these quantities sharply. 
In Microcanonical Thermodynamics we don't need von 
Neumann's definition (4) and can work on the level of 
Boltzmann-Planck's original definition of S by cqs.(l and 
3). We thus explore Statistical Mechanics and entropy on 
its most fundamental level. This has the great advantage 
that the axiomatic level is extremely simple. Because it 
does not demand scaling or cxtensivity, it can further be 
applied to the much wider group of non-extensive systems 
from nuclei to galaxies [4] and address the original object 
of Thermodynamics for which it was invented some 150 
years ago: phase separations. 

The Boltzmann-Planck formula has a direct and simple 
but deep physical interpretation: W or S are the mea- 
sure of our ignorance about the com,plete set of initial 
values for all 6N microscopic degrees of freedom which 
are needed to specify the N-body system unambiguously 
[5] . Usually we only know the few timc-indcpendcnt con- 
trol parameters E,N,V, - ■ ■ , of the system (conserved or 
very slowly varying) but have no control of the other fast 
changing degrees of freedom. If the information would be 
complete, i.e. all 6iV-degrees of freedom sharply known 
at some time to, W would be a single cell of size {2TTh)^^ 
in the SA^-dim. phase space, and S would be 0. 

However, many initial values are consistent with the 
same values of E,N,V,- ■ ■ , (redundancy of control pa- 
rameters). All the other degrees of freedom vary with 
time, usually fast, and are normally not under control. 
The manifold of all these points in the 6A^-dim. phase 
space is the microcanonical ensemble. It has a well 
defined geometrical size W and due to eq.(l) a non- 
vanishing entropy S{E, N,V, ■ ■ ■). The dependence of 
S{E, N,V,- ■ ■) on its arguments gives the complete ther- 
mostatics and equilibrium thermodynamics. 

It is further clear that the Hamiltonian (Liouvillean) 
dynamics of the system cannot create the missing in- 
formation about the initial values. I.e. the entropy 
S{E, N,V,- ■ ■) cannot decrease. As is further worked 
out in [6] and more recently in [4] the inherent finite res- 
olution of the macroscopic description implies usually an 
increase of or 5 with time when a constraint is re- 
moved (Second Law of Thermodynamics). 

It should be emphasized that this proper statistical Aci- 
inition of the entropy S{E,N,V,- ■ ■) characterizes the 
whole microcanonical ensemble [5] . It measures the total 
number of the possible states of the system under the in- 
formation given. Eqs. (2,3) are in clear contrast to the 



definition suggested by [7-10]: 

^^ui^i^) = J A.!(2.ft)3/ Q(^-^^'^'^» 

Eall cigcnstatcs n of H with given N,V, //^\ 
and a < < E 
= gSb^ik{E,N,V) (-7^ 

Our statistical definition (2,3) called Sgurj in [10] has a 
priori not much to do with the adiabatic invariants of a 
single trajectory in the 6A^-dim. phase space. 

As we will see later in section (III) the entropy S{E) 
for non-extensive systems can be convex in contrast to 
the assumption made by Hertz [7]. Then the tempera- 
ture or the mean kinetic energy per particle do not con- 
trol the energy flow during the equilibration of two sys- 
tems in thermal contact as was the crucial argument of 
Hertz [7] to introduce S'bu;fe,(5-7). In contrast, the Boltz- 
mann entropy eq.(2) does still control the direction of the 
energy-flow under equilibration even in these somewhat 
astonishing and unconventional though ubiquitous situa- 
tions. This clearly underlines the role of the Boltzmann 
entropy (2) as the fundamental quantity and not the tem- 
perature for statistical mechanics. I am fully aware that 
this contradicts opposite statements given by [7, 8, 10]. 

Moreover, the definition of the bulk entropy 
Sbuik{E, N,V) (5-7) contains also, energetically, inacces- 
sible states with less energy, which arc clearly excluded 
by our knowledge of the total energy. Thus, S^uik does 
not represent the redundancy, or ignorance, of the avail- 
able information. Therefore it is not acceptable by con- 
ceptual reasons. This is a clear warning that one should 
not follow some formal definitions and forget the original 
clear meaning of entropy and statistics as is often done 
in the literature. Of course, for macroscopic systems the 
difference between the two alternative definitions of the 
entropy eqs. (2-3) and eqs. (5,7) disappears. 

II. THE ZERO'TH LAW IN CONVENTIONAL 
EXTENSIVE THERMODYNAMICS 

In conventional (extensive) thermodynamics thermal 
equilibrium of two systems (1 & 2) is established by 
bringing them into thermal contact which allows free en- 
ergy exchange[27]. Equilibrium is established when the 
total entropy 

Stotai{E, E,) = Si{E,) + S2{E - El) (8) 

is maximal: 

dStotai{E,Ei)\E = dSi{Ei) + dS2{E-Ei)=0. (9) 

Under an energy fiux A£'2^i from 2 — > 1 the total en- 
tropy changes to lowest order in AE by 

AStotallE = iPl - f32)AE2^l (10) 

/3 = dS/dE = 1 (11) 
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Consequently, a maximum of Stotal{E = Ei + E2, Ei)\e 
will be approached when 

sign{AStotai) = sign(T2 - Ti)sign(A£;2^i) > (12) 

From here Clausius' first formulation of the Second Law 
follows: " Heat always flows from hot to cold" . Essential 
for this conclusion is the additivity of S under the split 
(eq.8). There are no correlations, which are destroyed 
when an extensive system is split. Temperature is an 
appropriate control parameter for extensive systems. 



III. NO PHASE SEPARATION WITHOUT A 
CONVEX, NON-EXTENSIVE S{E) 

Since the beginning of Thermodynamics in the middle 
of the 19. century its main motivation was the description 
of steam engines and the liquid to gas transition of water. 
Here water prefers to become inhomogcneoiis and develop 
a separation of the gas phase from the liquid, i.e. water 
boils. As conventional canonical statistics works only for 
homogeneous, infinite systems, phase separation remains 
outside of standard Boltzmann-Gibbs thermo-statistics. 

The weight e'^^^^"^/"^ of the configurations with en- 
ergy E in the definition of the canonical partition sum 

POO 

Z{T) = / e^(^)-^/^rf^ (13) 

becomes here bimodal, at the transition temperature it 
has two peaks, the liquid and the gas configurations 
which are separated in energy by the latent heat. Con- 
sequently S{E) must be convex and the weight in (13) 
has a minimum between the two pure phases. Of course, 
the minimum can only be seen in the microcanonical en- 
semble where the energy is controlled and its fluctua- 
tions forbidden. Otherwise, the system would fluctuate 
between the two pure phases by an, for macroscopic sys- 
tems even macroscopic, energy AE ^ Eiat of the order of 
the latent heat. I.e. the convexity of S{E) is the generic 
signal of a phase transition of first order and of phase- 
separation[l]. It is amusing that this fact that is essential 
for the original purpose of Thermodynamics to describe 
steam engines was never treated correctly in the past 150 
years. 

The ferromagnetic Potts-model illuminates in a most 
simple example the occurrence of a convex intruder in 
S{E) which induces a backbending caloric curve T{E) = 
{dS/dE)~^ with a decrease of the temperature T{E) with 
rising energy [11] [28]. A typical plot of s(e, N) = S{E = 
Ne)/N in the region of phase separation is shown in 
fig(l). Section IV discusses the general microscopic rea- 
sons for the convexity. 

This has far reaching consequences which are crucial 
for the fundamental understanding of thermo-statistics 
and Thermodynamics: Let us split the system of fig- 
ure (1) into two pieces a & 6 by a dividing surface, 
with half the number of particles each. The dividing 



FIG. 1: Ferromagnetic Potts model {q = 10) on a 50 * 50- 
lattice with periodic boundary conditions in the region of 

phase separation. At the energy ei per lattice point the sys- 
tem is in the pure ordered phase, at ea in the pure disordered 
phase. At Ca httle above ei the temperature Ta = 1//3 is 
higher than T2 and even more than Ti, at ei, a little below 
63. At Ca the system separates into a few bubbles of disor- 
dered phase embedded in the ordered phase or at ei, into a 
few droplets of ordered phase witiiin tiic disordered one. If we 
combine two equal systems: one with the energy per lattice 
site 60 = ei -I- Ae and at the temperature Ta, the other with 
the energy ej = — Ae and at the temperature Ti < Ta, 
and allowing for free energy exchange, then the systems will 
equilibrize at energy 62 with a rise of its entropy. The temper- 
ature Ta drops (cooling) and energy (heat) flows (on average) 
from b a. I.e.; Heat flows from cold to hot! Thus, the 
Clausius formulation of the Second Law is violated. This is 
well known for self-gravitating systems. However, this is not 
a peculiarity of only gravitating systems! It is the generic sit- 
uation at phase separations within classical thermodynamics 
even for systems with short-range coupling and has nothing 
to do with long-range interactions. 



surface is purely geometrical. It exists only as long 

as the two pieces can be distinguished by their differ- 
ent energy/particle Ca and e^. Constraining the energy- 
difference 66 — Ca = Ae between the two, reduces the 
number of free, unconstrained degrees of freedom and re- 
duces the entropy by —2ASsurf-corr.- (Moreover, if the 
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effect of the new surface would also be to cut some bonds: 
before the split there were configurations with attractive 
interactions across the surface which are interrupted by 
the division, their energy shifts upwards outside the per- 
mitted band- width eo, and thrown out of the partition 
sum (3). I.e. the entropy will be further reduced by the 
split.) 

If the constraint on the difference et — Ba is fully relaxed 

and gf, — Ca can fluctuate freely at fixed 62 = (e^ + eb)/2, 
the dividing surface is assumed to have no further phys- 
ical effect on the system. 

For an extensive system [S{E,N) = Ns{e = E/N) = 
2S{E/2,N/2)]. One would argue as follows: The com- 
bination of two pieces of N/2 particles each, one at 
Ca = 62 — Ae/2 and a second at e& = e2 + Ae/2, must lead 
to S{E2,N) > S{Ea/2, N/2) + S{Eb/2, N/2), the simple 
algebraic sum of the individual entropies because by com- 
bining the two pieces one normally looses information. 
This, however, is equal to [S{Ea, N) + S{Eb, N)]/2, thus 
S{E2,N) > [S{Ea,N) + S{Eh,N)]/2. I.e. the entropy 
S{E, N) of an extensive system is necessarily concave. 

For a non- extensive system we have in general 
S{E,N) > 2S{E/2,N/2) because again two separated, 
closed pieces have more information than their unifica- 
tion. Now, if E2 is the point of maximum positive cur- 
vature of S{E,N) (convexity = upwards concave like 
y = x^) we have S{E2,N) < [S{Ea,N) + S{Eb,N)]/2 
like in fig.(l). However, the r.h.s. is larger than 
S{Ea/2, N/2) + S{Eb/2, N/2). I.e. even though S{E, N) 
is convex at const. N, the unification of the pieces with 
Ea/2,N/2 and Eb/2,N/2 can still lead to a larger en- 
tropy S{E2,N). 

The difference between [S{Ea,N) + S{Eb,N)]/2 
and S{Ea/2, N/2) + S{Eb/2, N/2) wc call henceforth 
ASsurf-corr- The correct entropy balance, before and 
after establishing the energetic split eb > Co of the sys- 
tem, is 

Safter — Sbefore = ^ ASgurf-corr. — S2 < (14) 

even though the difference of the first and the last term 
is positive. 

In the inverse direction: By relaxing the constraint and 
allowing, on average, for an energy-flux [AEb^a > 0) op- 
posite to Ta — Tb > 0, against the temperature- gradient 



(slope), but in the direction of the energy-slope, the en- 
tropy Stotai ~^ S2 increases. This is consistent with the 
naive picture of an energy equilibration. Thus Clausius ' 
"energy flows always from hot to cold", i.e. the domi- 
nant control-role of the temperature in thermo- statistics 
[7] is violated. Of c;ourse this shows again that unlike 
to extensive thermodynamics the temperature is not the 
appropriate control parameter in non- extensive systems. 

In the thermodynamic limit N — > 00 of a sys- 
tem with short-range coupling ASsurf-corr. ~ N^^^, 

ASsurf-corr./N = ASsurf-corr. OC A/'"^/^ mUSt gO tO 

due to van Hove's theorem. 
IV. THE MICROCANONICAL MECHANISM 
LEADING TO CONDENSATION, PHASE 
SEPARATION AND THE ORIGIN OF THE 
CONVEXITIES OF S{E) 

Many theoretical applications of microcanonical ther- 
modynamics to realistic examples of hot nuclei, atomic 
clusters, and rotating astro-physical systems have been 
been presented during the past 20 years and show con- 
vex intruders in the microcanonical entropy and, conse- 
quently, negative heat capacities, c.f. e.g. the publication 
list on my WEB-site http://www.hmi.de/people/gross/, 
see also refs.[12 15]. Here I will illuminate the general 
microscopic mechanism leading to the appearance of a 
convex intruder in S{E, N,V,- ■ ■). 

I assume the system is classical and obeys the Hamil- 
tonian: 

N 2 
i 

i<j 

In this case the system is controlled by energy and vol- 
ume. 



A. Liquid— gas transition 

The microcanonical sum of states or partition sum is: 
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NT{m/2){27rh?y^/^ J^N [ E-Eo ) ^ ' 



Wid-gas {E-Eo,N,V)x Wint {E-Eo,N,V) 

f,[Sid-gas+Sint] 



(17) 



yN^^^(3N-2)/2 JjN ^3/2 

Wid-sas{E,N,V) = ^,r(3^/2)(2^;,2')3iv;2 (18) 

/ ^ \ (3iV-2)/2 

M^,„,(i;-iJo,iv,v) = y^^^e(£;-$»*[{r}])fi gj^^ M (i9) 



V is the spatial volume. Eq — min [{ r }] is the energy 
of the ground-state of the system. 

The separation of W{E, N, V) into Wid-gas and Wint 
is the microcanonical analogue of the split of the canoni- 
cal partition sum into a kinetic part and a configuration 
part, in the micro, however, without assuming any ther- 
modynamic limit or homogeneity: 



3N/2 



/ 



yJV 



(20) 



In the thermodynamic limit the order parameter of 
the (homogeneous) liquid-gas transition is the density. 
The transition is linked to a condensation of the system 
towards a larger density controlled by pressure. For a 
finite system we expect the analogue. However, here it 
is controlled by the constant available system volume V. 
At low energies the N particles condensate into a droplet 
with much smaller volume Vojv, the system does not fill 
the volume V. A'^ — 1 internal coordinates arc limited to 
Vojv- Only the center of mass of the droplet can move 
freely in V (remember we did not fix the CM. in eq.(16) ). 
The system does not fill the 3iV-configuration space V^. 

1 /3 

Only a stripe with width VqJy in 3(A^ — 1) dimensions 
of the total 3A/^-dim space and with the width V^^^ in 
the remaining 3 dimensions of center of mass motion is 
populated. The system is non-homogeneous even though 
it is equilibrized and, at low energies, internally in the 
single liquid phase. It is not characterized by an intensive 
homogeneous density. In fact, Wint{E — Eo,N, V) can be 



written as: 



Wint{E-Eo,N,V) = 

tN def 



V{E, N) 



V 



N 



< 1 



(21) 



[ViE,N)Y 

(3JV-2)/2 



E-Eo J 
Sint{E-Eo,N,V) = Nln 



V{E,N) 
V 



(22) 
< (23) 



The first factor e(£;-$™*[{r'}]) in eq(22) eliminates the 
energetically forbidden regions. Only the potential holes 
(clusters) in the 3A''-dim potential surface $*"*[{r}] re- 
main. Their volume V^{E,N) < is the accessible 
part of the 3A/'-dim-spatial volume where $'"*[{r}] < E. 
I.e. V'^(E,N) is the total 3A^-dim. eigen- volume of the 
condensate (droplets), with N particles at the given en- 
ergy, summed over all possible partitions, clusterings, in 
3A'^-configuration space. The relative volume fraction of 
each partition compared with {E, N) gives its relative 
probability. V^{E,N) has the limiting values: 



[V{E,N)f = <, „ ;v-l 



for E in the gas phase 
V for E = Eo 



W,nL {E - Eq, N, V) and Sint{E - Eq, N, V) have the lim- 
iting values: 

Wint{E - Eo) < 1, ^ Sint{E - Eo, N) < 
J 1 > 

^ U^]^"-^^ E^Eo 



(24) 



f E:^ 



(25) 
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All physical details arc encrypted in Wint{E—EQ, N, V) 
or Sint{E - Eo,N,V) alias Nln[V{E,N)], c.f. eqs.(21- 
25): If the energy is high the detailed structure of 
r*}] is unimportant WtnL ~ 1, StnL ~ 0. The sys- 
tem behaves like an ideal gas and fills the volume V. At 
sufficiently low energies only the minimum of [{ r }] 
is explored by WinL{E — Eo,N,V). The system is in a 
condensed phase, a single liquid drop, which moves freely 
inside the empty larger volume V, the 3{N — 1) internal 
dofs are trapped inside the reduced volume Vqn ^ V. 

One can guess the general form of N ln[V(-E, N)]: Near 
the groundstate E>Eq it must be flat « (A''— 1) ln[Vojv] + 
ln[V — Vojv] because the liquid drop has some eigen- 
volume Vojv in which each particle can move (liquid). 
With rising energy ln[V(i?, A^)] rises up to the point 
{Etrans) where it is possible that the drop fissions into 
two. Here an additional new configuration opens in 3A''- 
dim configuration space: Either one particle evaporates 
from the cluster and explores the external volume V, or 
the droplet fissions into two droplets and the two CM 
coordinates explore the larger V. This gives a sudden 
jump in Sint{E) by something like ~ ln{ ^ 

similar in the second case. 

Of course, this sudden opening of additional parts of 
the 3A^-dim. configuration space gives also a similar jump 
upwards in the total entropy 



V"^"""' } and 



S{E) 



Sid—qas ^" 



oc In 



(3iV-2)/2 



(26) 



by ~ ln{ 



V~Vo(]v-i) 
Vo(iv-i) 



}. Later further such "jumps" may 

follow. Each of these "jumps" induce a convex upwards 
bending of the total entropy S{E) (eq.l7). Each is con- 
nected to a bifurcation and bimodality of e^(^)~^/^ and 
the phenomenon of phase- separation. 

In the conventional canonical picture for a large num- 
ber of particles this is forbidden and hidden behind the 
familiar Yang-Lee singularity of the liquid to gas phase 
transition. 

In the microcanonical ensemble this is analogue to the 
phenomenon of multifragmentation in nuclear systems 
[1, 16]. This, in contrast to the mathematical Yang-Lee 
theorem, physical microscopic explanation of the liquid 
to gas phase transition sheds sharp light on the phys- 
ical origin of the transition, the sudden change in the 
inhomogeneous population of the SA'-dim. configuration 
space. 



B. Solid— liquid transition 

In contrast to the liquid phase, in the crystal phase 
a molecule can only move locally within its lattice cage 
of the size d^ instead of the whole volume Vqn of the 
condensate. I.e. in eq.(25) instead we have Sint 



C. Summary of section IV 

The gas-liquid transition is linked to the transition 

from uniform filling of the container volume V by the 
gas to the smaller eigen-volume of the system Vo in its 
condensed phase where the system is non-homogeneous 
(some liquid drops inside the larger empty volume V). 
First 3(A^— 1), later at higher energies less and less posi- 
tional degrees of freedom condensate into the drop. First 
three, then more and more dofs. (CM. -coordinates of 
the drops) explore the larger container volume V leading 
to upwards jumps (convexities) of Si,u{E). The volmne 
of the container controls how close one is to the criti- 
cal end-point of the transition, where phase-separation 
disappears. Towards the critical end-point, i.e. with 
smaller V the jumps ln[V — Vo] — ln[Vo] become smaller 
and smaller. In the case of the solid-liquid transition, 
however, the external volume V of the container confines 
only the CM. -position of the crystal resp. the droplet. 
As long as both d^ and Vqat remain small compared to V, 
the latter has no direct influence on Sint (E). The entropy 
jumps during melting are by ASini oc In Vqtv — Inrf^. 

At the surface of a drop > Ea = min*'"*, i.e. the 
surface gives a negative contribution to Sint in eq.(22) 
and to S at energies E>Eq, as was similarly assumed in 
section (III) and explicitly in eq.(14). 



V. ACKNOWLEDGEMENT 

I am grateful to J. Moller for insistent, therefore 
helpful, discussions. J.F. Kenney made numerous 
illuminating suggestions, I hope he finally can agree with 
this version. 



VI. APPENDIX, A SIMPLE MODEL FOR 
CONDENSATION 

Let's assume the various potential-pockets are attrac- 
tive square-wells with depths ^\ < 0. With the (some- 
what schematic) abbreviation: 



7(K) := 



j d^^r @{E - $'"*[{r }]) {e - $'"*[{r }]) 



(27) 

(3Af-K)/2 



iV-A 



A=l,Ar 



= $:((£;- $,)(3-'-K)/2\ V - ^ Vo. n 



k=l 



k=l 



Here $a < are different topological regions of ^>"* [{ r 
}]: Single cluster (A = 1), two clusters (A = 2), several 
clusters or, finally, if energetically possible, (A ~ N) free 
particles, $ = 0. Depending on A there are 3A CM. 
coordinates which can move over the whole volume V in 
contrast to the rest 3(A^ — A) internal coordinates which 
are limited to smaller cluster volumes Vok- 
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All particles condensed into one single cluster (A = 1) 
will have the lowest potential $i and consequently the 
largest — $. On the other hand the volume fac- 
tor ^ V/Vo(7v-i) ^ 1 appears only once. Then other 
terms with A > 1 may dominate if (V/Vo n-x)^ * {E — 

$^)(3JV-K)/2 > V/Vo(^-i) * {E - $,)(3iV-W)/2_ 

In terms of /(K) the partition sum W{E) and its 
derivatives are expressed as: 

W{E,N,V) oc 7(2) (28) 
aHV(E.N.V) ^ (3A>-2)(3iV-4)^^^^ 



dE^ 



3N -2 7(4) 
3N-2 



P'{E) 



3iV-4/(6) 3A^-2 fI{Ay 



1(2) 



To get a convexity (/?' > 0), we must have: 



mm) 



> 1+ 



3/^-4 



(31) 



(32) 



(33) 



2. Bimodal distribution, phase-separation 

When the massdispersion (27) is bimodal, when e.g. 
two terms in the sum dominate, i.e. when 



7(K) 



(38) 



(N-\ \ ^ N-X 
fc=i / fe=i 

then the different 7(K) arc not anymore proportional for 
different K as in eq.(34), and the curvature /?', eq.(32) 
can very well be positive. 

I.e. as anticipated in our previous discussion in section 
(IV), the volume- factor V/Vo > 1 controls whether the 
system can perform a phase-transition and whether the 
transition has a critical end-point and if, where it is, i.e. 
where the region of convexity disappears. 



3. Hard disks, the Alder- Wainwright Transition [17, 18] 



For hard disks the potential is: 
d -- 



m 

Eq 



r — r \ 
00 d< ddisk 
d > ddisk 



i>k 





rk) 



(39) 
(40) 
(41) 

(42) 



A. Some examples 



7(K) = 7;(3^-^)/2^(V-Voiv)X"^^^ 



(43) 



1. A single phase 

When calculating eq. (31-33) one immediately sees that 
for a sharply mono- dispersed sum (27), A ~ Aq (single 
phase) 



7(K-2) ~ {E-<^^,)I{K) 



(34) 



Then l.h.s. of criterion (33) becomes 1 and there is no 
convexity. 



I.e. the heat capacity is: 



3A^- 



1 



3iV- 2 



E-^Xo 



^ ^ 13' 2 



(35) 
(36) 

(37) 



equal to the ideal gas value, independently of whether 

the system is condensed or not. (Remember, in eq.(27) 
we assumed for simplicity that $a has a flat bottom). 



With Vojv = Nd^/V2 the closed packing volume. 

Due to the criterion (33) there is no bimodality vs. 
energy and the heat capacity is everywhere the ideal-gas 
value NciE) = 

The only interesting change of W is as function of the 
external volume V. To get some feeling for this, we make 
the following simplified ansatz: 

We assume each cluster with n disks has the closed- 
packing volume Von = -^Von- Moreover, each cluster 
has a constant surface layer of volume d around to sepa- 
rate it from other clusters. In a configuration of A clusters 
there are 3A positional degrees of freedom that can move 
free in volmiie of V — Vo — Ad whereas 3(A^ — A) internal 
positional degrees of freedom are limited to the conden- 
sate volume Vq^"'^^ For number of states we have after 
these simplifying assumptions, cf. eq.(43): 



1^ oc ^H^A 

A 

= (V - Von - Ad)^ x (Vow) 



N-X 



(44) 



The form of W\ as function of the external volume is 
shown in fig(2): 

Evidently, WaC^^) is for every V a sharp mono- 
dispersed function of number of clusters and shows no 
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bimodality v.s. A. Consequently there is no phase- 
separation also in the P{V) curve. This is not a phase- 
transition of first order. 



FIG. 2: Simplified Alder- Wainwright transition, Normalized 
cluster- number (A) distributions Wa(V) as function of the ex- 
ternal volume V. 



[1] D.H.E. Gross. Microcanonical thermodynamics: Phase 
transitions in "Small" systems, volume 66 of Lecture 
Notes in Physics. World Scientific, Singapore, 2001. 

[2] Alfred Wehrl. General properties of entropy. 
Rev.Mod.Phys., 50:221-260, 1978. 

[3] J. von Neumann. Gott.Nachr., 273:20-25, 1927. 

[4] D.H.E. Gross. A new thermodynamics from nuclei to 
stars. Entropy, 6:158 179. (2004). 

[5] J.E. Kilpatrick. Classical tlicnnostatistics. In H. Eyring, 
editor. Statistical Mechanics, number II, chapter 1, pages 
1-52. Academic Press, New York, 1967. 

[6] D.H.E. Gross. Ensemble probabilistic equilibrium and 
non-equilibrium thermodynamics without the thermo- 
dynamic limit. In Andrei Khrennikov, editor. Founda- 
tions of Probability and Physics, number XIII in PQ-QP: 
Quantum Probability, White Noise Analysis, pages 131- 
146, Boston, October 2001. ACM, World Scientific. 

[7] P. Hertz. Uber die mechanische Begriindung der Ther- 
modynamik II. Ann. Phys. (Leipzig), 33:537, 1910. 

[8] Victor L. Berdichevski. Thermodynamics of chaos and 
order. Longman, Edinburgh Gate, Harlow, England, 
1997. 

[9] Arnold Miinster. Statistical Thermodynamics. Springer, 
1969. 

[10] H.H. Rugh. Microthermodynamic formalism. 

Phys. Rev. E, 64:055101 1-4, 2001. 

[11] D.H.E. Gross, A. Ecker, and X.Z. Zhang. Microcanonical 
thermodynamics of first order phase transitions studied 
in the potts model. Ann. Physik, 5:446-452, 1996,and 
http:/ /xxx. lanl.gov/abs/cond- mat/9607150. 

[12] Ph. Cliomaz and F. Gulminclli. Energy correlations as 
thermodynamical signal for phase transitions in finite 
systems. Nucl. Phys. A, 647:153-171, 1999. 

[13] Ph. Chomaz, V.Dufiot, and F. Gulminelli. Caloric curves 
and energy fluctuations in the microcanonical liquid-gas 
phase transition, preprint, pages 1-4, 2000. 

[14] Ph.Chomaz, F.Gulminelh, and V.Duflot. Topology 
of event distribution as a generalized definition of 



phase transitions in finite systems. preprint, cond- 
mat/0010365:l-4, 2000. 

[15] Ph.Chomaz and F. Gulminelli. First order phase transi- 
tions: equivalence between bimodalities and the Yang- 
Lee theorem, preprint, cond-mat/0210456:l-10, 2002. 

[16] D.H.E. Gross. Microcanonical thermodynamics and sta- 
tistical fragmentation of dissipative systems — the topo- 
logical structure of the n-body phase space. Physics Re- 
ports, 279:119-202, 1997. 

[17] B. J. Alder and T.E. Wainwright. Phase transition in 
elastic disks. Phys. Rev., 127:359-361, 1962. 

[18] J.F. Kenney. The evolution of multicomponent systems 
at high pressures: I. the high-pressure, supercritical, gas- 
liquid phase transition. Fluid Phase Equilibria, 148:21- 
147, 1998. 

[19] L. G. Moretto, J. B. Elliot, L. Phair, , and G. J. Woz- 
niak. Negative heat capacities and first order phase 
transitions in nuclei. Phys. Rev.C, 66:041601 (R),nucl- 
thcor/0208024, 2002. 

[20] D.H.E. Gross. Statistical decay of very hot nuclei, the 
production of large clusters. Rep. Progr. Phys., 53:605- 
658, 1990. 

[21] M. D'Agostino, F. Gulminelli, Ph. Chomaz, M. Bruno, 
F. Cannata, R. Bougault, F. Gramegna, I. fori, 
N. le Neindre, G.V. Margagliotti, A. Moroni, and G. Van- 
nini. Negative heat capacity in the critical region of 
nuclear fragmentation: an experimental evidence of the 
liquid-gas phase transition. Phys.Lett.B, 473:219-225, 
2000. 

[22] B. Borderie. Dynamics and thermodynamics of the 
liquid-gas phase transition in hot nuclei studied with the 
indra array J. Phys. G: Nucl.Part.Phys., 28:R212-R227, 
2002. 

[23] M. Schmidt, R. Kusche, T. Hippler, J. Donges, W. Ko- 
rnmiiller, B. von Issendorff, and H. Haberland. Neg- 
ative heat capacity for a cluster of 147 sodium stoms. 
Phys. Rev. Lett, 86:1191-1194, 2001. 

[24] M. Schmidt and H. Haberland. Caloric curve across the 



9 



liquid-to-gas change for sodium clusters. Pys. Rev. Lett., 
87:203402- 1 4, 2001. 

[25] D.H.E. Gross. Comment on "negative heat and first or- 
der phase transitions in nuclei" by Moretto et al. pages 
http://arXiv.org/abs/nucl-th/0208034, 2002). 

[26] F. Gulminelli and Ph. Chomaa. Comment on "negative 
heat capacities and first order phase transitions in nuclei" 
by l.g. Moretto et al. Caen preprint, LPCC 99-17, pages 
nucl-th/0304058, 2003. 

[27] This and the next (III.) section discuss mainly systems 
that have no other macroscopic (extensive) control par- 
rameter besides the energy. E.g. the particle density is 



not changed and there arc no chemical reactions. 
[28] This backbending of T{e) has nothing to do with the 
periodic boundary condition used here, as claimed erro- 
neously by Moretto et.al [19]. Many realistic microcanon- 
ical calculations for nuclear fragmentation as also for 
atomic-cluster fragmentation, c.f. [1, 20], show the back- 
bending without using periodic boundaries. These are ver- 
ified by numerous experimental data in nuclear physics 
c.f. [21, 22] and cluster physics [23, 24]. The errors of the 
above paper by Moretto et. al. are commented in more 
detail also in [25, 26] 



A A 

:V-VO-Kd) /VO 



V0 = 4.6: 



d = 0.03 



V=5 V=1 V = 1 5 



t 

t 
t 
t 



I 

T 

t 

T 
i 

T 

i 
t 



1 



♦ 

♦ 
♦ 
♦ 
♦ 
♦ 
♦ 
♦ 
♦ 
♦ 
♦ 
♦ 
♦ 



♦ 
♦ 
♦ 
♦ 
♦ 
♦ 
♦ 
♦ 
♦ 
♦ 
♦ 
♦ 
♦ 



V=20 
ft 



an 



V=30 



TT 

f T 

▼ T 

▼ T 



T 

t 
t 
f 
t 
t 

T 
T 

T 

▼ 
T 

t 

T 



t 

T 

T 

t 

t 

▼ 

t 

t 

T 

T 

t 
T 
T 
t 



LA 



100 



300 



400 



500 



A 



